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ABSTRACT 

Infrared Dark Clouds (IRDCs) are ideal laboratories to study the initial processes of high-mass star and star cluster formation. We 
investigated star formation activity of an unexplored filamentary dark cloud (size ~ 5.7 pc x 1.9 pc), which itself is part of a large 
filament (~ 20 pc) located in the S254-S258 OB complex at a distance of 2.5 kpc. Using Multi-band Imaging Photometer (MIPS) 
Spitzer 24 pm data, we uncover 49 sources with signal-to-noise ratio greater than 5. We identified 45 sources as candidate young 
stellar objects (YSOs) of Class I, Flat-spectrum, and Class II nature. Additional 17 candidate YSOs (9 Class I and 8 Class II) are 
also identified using JHK and Wide-field Infrared Survey Explorer (WISE) photometry. We find that the protostar to Class II sources 
ratio (~ 2) and the protostar fraction (~ 70 %) of the region are high. When the protostar fraction compared to other young clusters, 
it suggests that the star formation in the dark cloud was possibly started only 1 Myr ago. Combining the near-infrared photometry of 
the YSO candidates with the theoretical evolutionary models, we infer that most of the candidate YSOs formed in the dark cloud are 
low-mass (< 2 M©) in nature. We examine the spatial distribution of the YSOs and find that majority of them are linearly aligned along 
the highest column density line (N(H 2 )~ 1 X 10^^ cm"^) of the dark cloud along its long axis at mean nearest neighbor separation 
of ~ 0.2 pc. Using observed properties of the YSOs, physical conditions of the cloud and a simple cylindrical model, we explore 
the possible star formation process of this filamentary dark cloud and suggest that gravitational fragmentation within the filament 
should have played a dominant role in the formation of the YSOs. From the total mass of the YSOs, gaseous mass associated with 
the dark cloud, and surrounding environment, we infer that the region is presently forming stars at an efficiency ~ 3% and a rate ~ 30 
Mq Myr"^ and may emerge to a richer cluster. 

Key words, infrared, young stellar objects, dark cloud, filament, star formation; individual: S254-S258 complex 


1. Introduction 

Recently, filaments and infrared dark clouds (IRDCs) have re 
ceived special attention because they are the potential progen 
itors of cluster formation (e.g., Rathborne et al.||2006[ Beuther 

& Sridharan 2007trPeretto & Fuller|2009t|Myers|2009[ [Henning 

et al.|20I0[ ' Battersby et al.|20II[|Ragan k al.|20I2t|Longmore 


et al.|20I27|Russeil et al.|20I3| ). Results from the Herschel S pace 

Telescop e, particularly the Gou ld Belt ( [Andre et al.||20I0| ) and 
Hi-GAL ( jMolinari et al.|[MT0| Legacy Surveys, have also em 
phasized the crucial role of filamentary structures on star and 
cluster formation. However, due to the complex nature of the 
interstellar medium (ISM), the modes of fragmentation and the 
physical processes that govern star formation in filaments are 
still under active discussion (e.g.,|Heitsch et al.|2009t[Wzquez 


Semadeni et al.|20Q9t|Pon et al.|20IIl|Andre|20I3t|Smith et al 

20I4| ). The identification and characterization of young proto¬ 
stars represent a key point to partly answer some of these ques¬ 


tions, as they still carry the imprint of the fragmentation of the 
primordial cloud. However, filaments and IRDCs are character¬ 
ized by high column densities and, therefore, the identification 
of such sources is often difficult. In the last decade, observa¬ 
tions at longer infrared wavelengths with Mid-infrared Imaging 
Photometer (MIPS; [Carey et aL 2009[) and Wide-field Infrared 
Survey Explorer (WISE; [Wright et al.|20I0[ ) allowed us to over¬ 
come the high extinction in such regions and explore the proper¬ 
ties and distribution of associated young sources. In this context, 
we present an investigation of a filamentary dark cloud located 
in the S254-S258 OB complex. 

The S254-S258 complex, situated a t a distance of 2.5 kpc, 
is part of the Gemini OB association ( [Chavarria et al.[[2008t 
[Bieging et al.[[2009l [Ojha et al.[[20II|). The complex contains 
six HII regions (S254, S255, S256, S257, S258 and S255B) pro¬ 
jected on or in the close vicinity of a long (~ 20 pc) filamentary 
cloud (see Eig. I). The search for young stellar objects (YSOs) 
in the complex has been carried out by [Chavarria et al.[ ( [2008[ ). 
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et al.lpOlT] ), and |Mucciarelli et aL] ( [2011| ), but these studies 
mainly focus on the central part of the complex, i.e., roughly over 
<13 square arcmin area around the massive young cluster S255- 


IR ( [Wang et al.||MTT| |Qjha et al.|[MTT| [Zinchenko 

located between the two evolved H ii regions S255 and S257. 
These studies have shown that the S254-S258 complex contains 
a rich population of Class I, Class II, Class III, and near-infrared 
(NIR)-excess sources formed in groups, and also in distributed 
mode across the complex. 

The present study focuses on the star formation activity of 
a part of the long filament (centered at 0^2000 = 06^13"^47^, 
^2000 = 17°54'37" or /=192°.76, /?=00°.10). The filamentary 
area under investigation is of size ~ 5.7 pc in length and 1.9 pc 
in width (marked with a rectangular box in Fig. 1). The region is 
highly obscured at optical wavelength regime, and corresponds 
to a dense part of the long filament. The region is dense as re¬ 
vealed by the presence of CS (J=2-l) emission, which is a tracer 
of dense (> 10"^ cm“^) gas. The distribution of CS emission (see 
Fig. 4 of [Carpenter et al.||1995| ) is also seen elongated in shape 
(size ~ 5.6 pc X 1.7 pc), with its long axis along the long axis 
of the large filament. IRDCs are dense molecular clouds seen as 
extinction features against the bright background. Molecular line 
and dust continuum studies of IRDCs have shown that they are 


cold (T < 25 K), dense (A^(H 2 ) ^ 10^^ cm ^), and massive (- 


10^-10^ Mq) structures with sizes 1-15 pc ([Carey et al. 


1998 


Simon et al.|2006[ ). Like other IRDCs (e.g., [Simon et al.|2006[ ). 


the area under investigation is dark in optical (e.g., in DSS2 Sur¬ 
vey images) as well as in mid-infrared (e.g., in MSX Survey im¬ 
ages; |Pnc££r^(|20^), and also cold (T < 15 K) and dense 
(A^(H 2 ) ~ 10^^ cm ^) in nature (discussed in Sects. 3.5 and 4.2). 


Hence, following the general nomenclature of the IRDCs ( Rath-[ 
borne et al.[2006[[Peretto & Fuller[200^ , hereafter, we designate 
the dense cloud under investigation as “IRDC G192.76-r00.10”. 

Massive OB2 stars ionize the surrounding ISM and create 
H II region around them. During expansion, H II region drives 
the shock front preceding ionization front. The shock front 
sweeps up the ISM to form a cold dense shell and star forma¬ 
tion can be induced by the instability of the shell ( e.g., “col¬ 
lect and collapse“ sc enario; [Elmegreen & Lad^[1977[). In the 
S254-S 258 complex, [Qjha et al.[ ( [2011[ ) and [Mucciarelli et al 


( 2011[ ) suggested the presence of induced second generation 


star formation at the peripheries of S255 and S257. Similar to 
here, evidences of second generation triggered star formations 
have been observed at the peripheries of several H ii regions 


(e.g., 

Zavagno et al.|2007l[Indebetouw et al.|2007[ 

Chu & Cm- 

endl[ 

2008t IChen et al. |2011l |Brand et all 2011 

Deharveng 

et al. 

[20121 [Getman et al.[[2012t [Jose et al.[[2013l 

Samal et al. 


2014 ), suggesting that the presence of an H ii region can infiu- 


ence star formation processes of a complex. However, in the 
S254-S258 complex, the IRDC G192.76-r00.10 region is rel¬ 
atively isolated from the feedback effects, as it is situated far¬ 
ther from the known evolved H ii regions (e.g., S254, S255, and 
S257). Moreover, shell-like filamentary structures resulted due 
to sweeping and compression of expanding H ii regions are gen- 
erally located parallel to the H ii region’s ionization front (e.g., 
[Deharveng et al.[[201Q[ [Jose et aL][2013[ ), whereas the filament 
axis of the IRDC G192.76-r00.10 region is perpendicular to 
the ionization front of the evolved H ii regions of the complex. 
Also the IRDC G192.76-h00.10 area does not seem to be con¬ 
taminated by polycyclic aromatic hydrocarbon (PAH) emissions 
at 11.3 yum and 12.7 pm of WISE 12 pm band (not shown); 
PAH emissions are generally found at the interaction zone of 
HII region and molecular cloud (e.g., [Samal et al.|2007[[Deh^ 


veng etal. 2010). These evidences suggest that the star formation 
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in the IRDC G192.76-1-00.10 region is unlikely to be influenced 
by the evolved H ii regions of the complex. 

Star formation in the IRDC G192.76-1-00.10 region is un¬ 
explored so far; only partly covered by [Chavarria et al.[ ( [2008[ ) 
at Spitzer-IRAC wavelengths (i.e., 3.6-8.0 pm) in their YSOs 
search. In the present work, along with data, we make 

use of MIPS, WISE, and NIR-JHK band data sets to study the 
star formation activity of IRDC G192.76-h00.10. We present our 
work with the following layout. We describe the observations 
and data reduction techniques in Sect. 2. In Sect. 3, the observa¬ 
tional results are presented, which include morphology of the 
region and identification, distribution, and characterization of 
YSOs. We discuss the properties and star formation processes 
of the IRDC G192.76-h00.10 cloud and its overall relation with 
the main filament in Sect. 4. 


Chavarria et al.[([2008 


2. Observations and Data Reduction 


The 24 pm observations of the IRDC G192.76-i-00.10 re¬ 
gion were downloaded from the Spitzer archive (Program ID 
20635), and cover an area ~ 7'0 x 8'0, centered at a 2 ooo = 
06'’13'"47", (52000 = 17°53'40". We downloaded the corrected 
basic-calibrated data (CBCD) images and the corresponding un¬ 
certainty files and performed source detection and photometry 
using the MOPEX/APEX software. To extract the flux, we ap¬ 
plied the APEX point-response function (PRF) fitting method to 
the detected sources. We used the zero-point value of 7.17 Jy 
from the MIPS Data Handbook to convert flux densities to mag¬ 
nitudes. We considered only those 49 sources whose signal-to- 
noise (S/N) ratio was found to be greater than 5. The magnitudes 
of these 49 sources along with their corresponding uncertainties 
are given in Table 1. We point out that the reported uncertainties 
are lower limit to the actual values, as the uncertainty in the ab¬ 
solute flux densities at 24 pm is ~ 4% ( [Engelbracht et al.|2007[ ). 

We also downloaded MIPS 70 pm CBCD image from the 
Spitzer archive (Program ID 20635), and performed photometry 
using APEX-PRF fitting method as discussed above. The area 
covered by the MIPS 70 pm image is ~ 5'7 x 2'0, centered at 
0^2000 = 06M3"^47^ ^2000 = 17°54'45". We detected only four 
point sources with S/N > 5. Zero-point value of 0.778 Jy, adopted 
from the MIPS Data Handbook, was used to convert the PRF 
fitted flux densities of these sources to magnitudes. 

Deep NIR observations of the IRDC G192.76-1-00.10 region 
in / (T = 1.25 pm), H {A = 1.63 pm\ and K {A = 2.1 4 pm) 
bands were obtained from the WFCAM Science Ar chive ([Ham- 
bly et al.[2008[ ), taken with the WFCAM instrument ( [Casali et ST 
2007[ ) at the United Kingdom Infrared Telescope (UKIRT). We 
performed photometry on the retrieved stacked images produced 
by WFCAM pipeline at Cambridge Astronomical Survey Unit 
(CASU) for an area of ~ 15' x 12' around the IRDC region cen¬ 
tered at (12000 = 06n3"^48^ & ^2000 = 17°52'50". Photometry 
on the images was done using the PSF algorithm of DAOPHOT 
package (Stetson 1987) in IRAK The PSF was determined from 
the bright and isolated stars of the field. For photometric cali¬ 
bration, we used isolated Two Micron All Sky Survey (2MASS) 
point sources ( [Cutri et al.|2QQ3[ ) having error <0.1 mag and rd- 
fiag “123". Rd-flag values of ’ 1’, ’2’ or ’3’ generally indicate the 
best quality detections, photometry and astrometry, respectively. 
A mean calibration dispersion of < 0.07 mag is observed in each 
band, indicating that our photometry is reliable within ~ 0.07 
mag. Saturated sources in our catalog were replaced by 2MASS 
sources. 

Only 9 out of the 49 MIPS detected sources have a coun¬ 
terpart in the NIR and Spitzer-\RAC bands from the catalog of 
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Fig. 1. Colour composite image of the complex obtaine d using the (110 GHz) column density map in red (from Chavarria et al. ( |2008| )), 
radio emission in green (from NVSS survey at 1.4 GHz; [Condon et al~ ( |1998| )) and optical emission in blue (from DSS2 survey). The abscissa 
(RA) and ordinate (Dec) are in J2000 epoch. North is up and east is left. The H ii regions discussed in the text are marked in white circles. The 
rectangular box represents the area studied in this work (the corresponding image at 24 fim is shown in Fig. 2). 


Chavarria et al.| ( [200^ . We thus mostly used WISE and our NIR 


point source catalogs to classify YSOs and for the construction 
of their spectral energy distributions (SEDs). The WISE survey 
( [Cutri & et al.||2012| ) provides photometry at four wavelengths: 
3.4, 4.6, 12 and 22 yum, with an angular resolution of 6'.'1, 6'.'4, 
6'.'5 and 12''0, respectively. We merged the MIPS catalog with 
the WI SE and NIR catalogs using a matching radius of 3'.'0 (fol¬ 
lowing Koenig et~ar|(2012| )). 

[Chavarria et al.| ( |2008| ) conducted observations in the J=l- 0 
spectral lines of ^"^CO and ^^CO using the Eive College Radio 
Astronomy Observatory (ECRAO) 14 m telescope. The ECRAO 
beam size is 45" in ^^CO and 46" in ^^CO. We used the ^^CO 
column density map of [Chavarria et al.| ( |2008| ) to study the gas 
content of the region. 


3. Results 

3.1. Morphology 

The MIPS 24 pm image of the IRDC G192.76-h00.10 region 
is shown in Pig. 2. The image displays a significant number 
of point-like sources aligned roughly in a linear sequence from 


south-east to north-west, and most of the point sources seem to 
be bounded by an elongated structure of ^^CO gas. Since dark 
cloud and dense molecular gas are the sites of new star forma¬ 
tion, these sources are possibly young protostars in their early 
evolutionary stages. In order to study the star formation activ¬ 
ity of the region, it is necessary to characterize and discuss the 
nature of these sources. 


3.2. Identification of Young Stellar Objects 


The circumstellar dust emission from the disk and infalling enve¬ 
lope of young stars gradually disappears with time as a function 
of their evolutionary phases. A number of classification methods 
are employed in the literature to classify evolutionary phases of 
YSOs; often used are the near- to mid-infrared spectral index a 
( Lada|198'7] ), the ratio of submillimeter to bolometric luminosity 
Lsubmm/Lboi (^.g., Andre et al.|2000| ), and the bolometric temper¬ 
ature Tboi (e.g., [Myers et al.|1998| ). Generally, submillimeter lu¬ 
minosity is taken to be the integrated luminosity at wavelengths 
T > 350 pm. The bolometric luminosity is calculated by integrat¬ 
ing the SED over extensive wavelength coverage, covering far- 
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Fig. 2. SpitzerMlPS 24 //m image of IRDC G192.76+00.10. The image 
has a field of view ~7:0 x 8^0 centered at a 2 ooo = 06^13'”47% 62000 = 
17°54'37" or /=192°.76, Z?=00°.10. The column density contours 
are also shown (contour levels are at 0.7, 0.9, 1.1, 1.4 x 10^^ cm“^). 



Fig. 3. The [3.4] - [24] vs. [3.4] colour-magnitude diagram for the de¬ 
tected MIPS 24 pm sources. The dotted lines separate the regions of 
Class I, Flat-spectrum. Class II, and Class III objects. The dashed lines 
denote the dividing line between the region occupied by contaminated 
sourc es (galaxies and disk-less stars) and YSOs (see also |Rebull et al.| 
|2011| ). The sources for which we do not have WISE 3.4 pm detection 
are represented in red. For these sources, we consider the magnitude of 
the faintest 3.4 pm counterpart of our 24 pm detections as an upper- 
limit. The arrow represents the direction of their colours, thus they are 
likely to be Class I YSOs. 


Table 1. MIPS 24 pm photometry for the 49 sources detected in the 
IRDC G192.76+00.10 region with a S/N ratio greater than 5 


ID 

RA (deg) 
J2000 

DEC (deg) 
J2000 

[24] yum 
mag 

a 

1 

93.424041 

17.944051 

8.568 + 0.051 

-0.43 

2 

93.416244 

17.943596 

7.645 + 0.032 

-0.29 

3 

93.386224 

17.934850 

5.809 + 0.017 

0.43 

4 

93.483044 

17.933043 

8.339 + 0.044 

-3.06 

5 

93.433941 

17.926279 

8.407 + 0.034 

-0.82 

6 

93.446990 

17.921782 

8.060 + 0.024 

-0.02 

7 

93.442216 

17.921563 

6.423 + 0.016 

-0.18 

8 

93.436925 

17.918236 

4.256 + 0.011 

0.93 

9 

93.424662 

17.917668 

6.554 + 0.017 

-0.31 

10 

93.447529 

17.914416 

4.802 + 0.011 

1.83 

11 

93.432413 

17.912164 

8.354 + 0.032 

-1.17 

12 

93.460228 

17.911043 

3.570 + 0.011 

1.92 

13 

93.457060 

17.906121 

8.200 + 0.029 

— 

14 

93.466088 

17.904649 

7.684 + 0.018 

-0.39 

15 

93.453099 

17.903396 

5.898 + 0.014 

-0.34 

16 

93.461708 

17.903433 

5.564 + 0.013 

1.27 

17 

93.441803 

17.901178 

8.424 + 0.033 

-0.85 

18 

93.457880 

17.899842 

9.098 + 0.065 

-2.51 

19 

93.468669 

17.895977 

6.052 + 0.014 

-2.61 

20 

93.476514 

17.889089 

6.750 + 0.010 

-0.64 

21 

93.503164 

17.888787 

4.377 + 0.012 

1.37 

22 

93.494384 

17.886028 

7.768 + 0.028 

0.34 

23 

93.483423 

17.882087 

8.463 + 0.057 

-0.55 

24 

93.498257 

17.880100 

3.893 + 0.001 

1.10 

25 

93.500148 

17.878929 

5.609 + 0.006 

0.52 

26 

93.496984 

17.877505 

5.850 + 0.008 

0.26 

27 

93.459530 

17.877001 

8.918 + 0.073 

— 

28 

93.482464 

17.875690 

6.792 + 0.017 

-0.13 

29 

93.392728 

17.875741 

7.388 + 0.032 

-2.62 

30 

93.482421 

17.872994 

7.248 + 0.026 

— 

31 

93.478138 

17.863109 

8.028 + 0.050 

-0.29 

32 

93.490105 

17.861847 

6.729 + 0.016 

-0.46 

33 

93.491665 

17.856101 

5.204 + 0.004 

2.08 

34 

93.481174 

17.852611 

8.152 + 0.061 

-0.32 

35 

93.497793 

17.852376 

5.718 + 0.017 

— 

36 

93.496107 

17.850057 

5.564 + 0.016 

0.26 

37 

93.493245 

17.848029 

6.739 + 0.019 

— 

38 

93.479544 

17.894587 

9.160 + 0.091 

— 

39 

93.474166 

17.930534 

9.010 + 0.075 

— 

40 

93.467083 

17.908493 

9.617 + 0.131 

— 

41 

93.477897 

17.937458 

9.523 + 0.140 

— 

44 

93.492631 

17.906890 

4.921 + 0.012 

0.17 

45 

93.437399 

17.921868 

6.806 + 0.018 

-0.56 

46 

93.450619 

17.906530 

7.407 + 0.035 

— 

47 

93.471334 

17.903989 

6.356 + 0.176 

1.83 

48 

93.444723 

17.913772 

7.178 + 0.250 

0.36 

49 

93.450026 

17.905202 

6.620 + 0.180 

1.10 


infrared and millimeter domain. Similarly, the bolometric tem¬ 
perature is defined as the temperature of a black body with the 
same mean frequency as the source SED over wide wavelength 
coverage. Since, the identified point sources typically have flux 
measurements at T < 24 pm, accurate determinations of LtoP 
Thai, and Lsubmm is not feasible with the data; thus, in this work, 
we classifled the sources based on their a values. 


Out of the 49 MIPS detections (see Sect. 2), 38 have WISE 
counterparts in all the four WISE bands. We estimated spectral 
index (a = d \og(AF^)ld log(/I), where E^ is the flux as a func¬ 
tion of wavelength, A) for these 38 sources from a linear fit to the 
fluxes in the range 3.4 to 24 pm (four WISE bands and MIPS 24 
pm band), and then classifled sources as Class I (a > 0.3), Elat- 
spectrum (0.3 >a> -0.3), Class II (-0.3 > tr > - 1.6), and Class III 
(a < -1.6) YSOs following [Evans et ^ ( |2009| ). Using the above 
approach, we And that out of the 38 sources, 16, 7, 11, and 4 
sources are Class I, Elat-spectrum, Class II, and Class III YSOs, 
respectively. The a values of the 38 sources are also tabulated in 
Table 1. Here we would like to mention that, though a is one of 
the most commonly used methods for the classiflcation of YSOs, 
it is highly susceptible to disk geometry and source inclination 
(e.g., |Robitaille et al.|2006[|Crapsi et al.|2008a| ). Also in this ap¬ 
proach the distinction of the Class 0 YSOs from the Class I is not 
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possible as no well-defined a criteria exists for Class 0 sources. 
This is because, Class 0 spectrum has lowest fiux densities from 
1.25 - 24 jum as expected for deeply embedded sources with 
massive, extincting envelopes, thus were largely not identified in 
the mid-IR prior to Spitzer. However, it is possible to distinguish 
them from the Class I sources based on their Lsubmm/Lboi and T^oi 
estimations as discussed in | Andre et al.|(|20001) and|Myers et al.| 
( [T9981 ). 


It has been found that a strong correlation between the colour 
and spectral index exists, which provides an acceptable proxy to 
classify YSOs (see, e.g., |Rebull et al. [2007112010| ). We used the 
[3.4] vs. [3.4] - [24] colour-magn itude diagram (see Fig. 3) to 
classify YSOs ( |Guieu et~ar]|2Q10| ). In Fig. 3, the regions occu¬ 
pied by Class Ill/main-sequence (MS), Class II, Flat-spectrum, 
and Class I objects are indicated. In this diagram, 34 sources are 
located in the regions occupied by Class I, Flat-spectrum, and 
Class II YSOs, whereas 4 sources are found in the zone of disk¬ 
less stars. We confirmed the nature of these four sources as disk¬ 
less stars with SED modeling (see Sect. 3.6); these sources will 
not be considered as YSOs in the followings. For 11 MIPS de¬ 
tections, 3.4 pm counterparts are not available. For these sources 
we then considered the magnitude of the faintest 3.4 pm coun¬ 
terpart of the 24 pm detections as an upper limit, in order to 
determine their positions on the [3.4] vs. [3.4] - [24] diagram. 
The approximate positions of these 11 sources are shown (red 
dots) in Fig. 3; they appear to be Class I YSOs. Out of these 
11 sources, 4 sources have K-band detection. We find that they 
also fall in the Class I zone in the K vs. K-[24] diagram (e.g.. 
Rebull et al.|2007| ), suggesting that these sources are likely to be 


Class I YSOs. For 4 Class I YSOs, we ha ve 70 pm fiux, we th us 
used the [24] vs. [24]-[70] diagram (e.g., |Rebull et al.|20Q7| ) to 
be more precise about their nature; their positions on the [24] vs. 
[24]-[70] diagram suggest that they are indeed Class I YSOs. 

In Fig. 3, most of the sources are having [3.4]-[24] colour 
> 4. A normal MS star has [3.4]-[24] colour around 0.0. Thus, 
a normal MS star would require a foreground visual extinction 
(Av ) ~ 210 mag to redden up to the location [3.4]-[24] ~ 4, 
whereas the mean Ay of the region is ~ 10 mag (see Sect. 3.3), 
suggesting most of the MIPS sources are most-likely YSOs with 
disks and/or envelopes. However, the YSO sample can be con¬ 
taminated by background dusty active galactic nuclei (AGN) 
and star-forming galaxies, as they have colours similar to that 
of YSOs. Hence, in Fig. 3, we marked the approximate zones 
of such sources based o n the Spitzer Wide-are a Infrared Extra- 
galactic Survey catalog ( [Lonsdale et al.|[2003] ) observations of 
the ELAIS N1 extra-galactic field. Two sources of our sam¬ 
ple (i.e., sources with [3.4]-[24] colour between 6 and 7 mag) 
fall in the extra-galactic zone. However, these sources are found 
close to other YSO candidates, and are likely members. Similar 
to galaxies, the colours of the asymptotic giant branch (AGB) 
stars can also mimic with the colour of YSOs. AGB stars show 
a steep spectral index at long wavelengths. Based on their SEDs, 
[Robitaille et al.j ( |2008j ) used a criterion [8.0] - [24] < 2.5 to 
identify them. In the absence of IRAC 8.0 pm data, we are un¬ 
able to apply [Robitaille et al.j ( [2008[ ) criterion to eli minate such 
sources. However, based on AGB stars distribution, [Robitaill^ 
et al.[ ( [2008[ ) established an empirical relationship that roughly 


reflects the expected number of AGB stars per square degree area 
on the sky for a given Galactic longitude and latitude. Using [Ro^ 
bitaille etn^ ( [2008[ ) relationship, we calculated that the AGB star 


contamination to our YSO sample is likely to be less than one. 

Above discussions suggest that the contamination of non- 
YSO sources to our MIPS identified YSO sample should be neg¬ 
ligible. We therefore considered all the 45 MIPS identified YSO 


candidates of different classes (27 Class I, 7 Flat-spectrum, and 
11 Class II) for further analyses. 



([ 3 . 4 ]-[ 4 . 6 ]), 


Fig. 4. Intrinsic K - [3.4] vs. [3.4] - [4.6] colour-colour diagram show¬ 
ing the distribution of field stars/MS sources (cross symbols), Class II 
YSOs (blue stars) and Class I YSOs (red stars). Slanted line indicates 
the boundary between the Class I and Class II YSOs. The control field 
sources are marked with open circles. The error bars in the bottom right 
comer show average errors in the colours. 


3.3. Additional YSO candidates 


The YSO detection using the above method is primarily limited 
by the 24 pm detection limit. Thus, we may be missing faint 
YSOs of the region. To overcome this problem, we matched 
the WISE catalog to our ^-band point-source catalog, and then 
selected YSO candidates using the intrinsic K - [3.4] vs. [3.4] 
- [4.6] colour-colour diagram as suggested by [Koenig et al. 
( [2Q12[ ). The intrinsic K - [3.4] vs. [3.4] - [4.6] diagram is shown 
in Fig. 4, and the zones of Class II and Class I YSOs are also 
marked. To construct intrinsic colour-colour diagram, we mea¬ 
sured the visual extinction towards individual YSOs from the 
H 2 colu mn density map (see Sect . 3.5), and then used e xtinction 
laws of [Bohlin et al.[ ( p^78 ) and [Flaherty et al.[ ( [2007[ ) to com¬ 
pute their dereddened K - [3.4] and [3.4]-[4.6] colours. We then 
used [Koenig et al.[ ( [2012[ ) criterion to select the Class I and Class 
II YSOs from the intrinsic K - [3.4] vs. [3.4] - [4.6] diagram. 
This helped us to identify 17 additional YSO candidates, which 
include 8 Class II and 9 Class I YSOs. To quantify the contam¬ 
ination of other sources to this selection, we selected a control 
field devoid of CO gas and located outside the filamentary area, 
but within the field of view of our NIR observations. We then 
looked for the distribution of control field sources on K - [3.4] 
vs. [3.4] - [4.6] diagram. These sources are also shown in Fig. 4 
(marked with open circles). Their distributions suggest that the 
majority of them are field stars, indicating that the contamination 
of other sources to the YSOs sample, selected based on K -[3.4] 
vs. [3.4]-[4.6] diagram, should only be a few. 

In summary, we have identified 62 YSO candidates in the 
IRDC G192.76-1-00.10 region with excess IR emission. Of these, 
43 have excess consistent with Class I plus Flat-spectrum YSOs, 
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and 19 have excess consistent with Class II YSOs. We note 
that even though the near to mid-infrared colours are very use¬ 
ful to identify YSOs, the YSO selection and classification can 
be biased due to disk geometry and/or source inclination along 
the line of sight. Infrared spectroscopic observations provide 
many useful indicators that are helpful for the YSO confirma¬ 
tion and disentangleme nt of deeply embedded protostars from 
the evolved YSOs (e.g.,|Reach|2007[ Spezzi et al.|20Q8t[Crapsi| 
|et al.|[2008b{ |Connelley & Greene 2010| ), however foreground 
absorption and edge-on disks can still co nfuse the YSOs clas¬ 
sification (e.g., [Pontoppidan et al.]|2005| ). Nevertheless, it has 
been found that the YSO selection scheme based on photometric 
colours provides a good representation of YSOs identification 
in star-forming complexes. As an example, [Spezzi et ak] ( |2QQ8| ) 
with spectroscopic study of Chameleon II cloud fou nd that 96% 
of the YSOs identified in Spitzer c2d Legacy survey ( [Evans etak 


20031 using photometric colours were true members. 



Fig. 5. Completeness of the 24 pm band data for the IRDC 
G192.76+00.10 region from artificial stars experiment. 



14 16 18 20 

K 


Fig. 6. The cumulative distribution of all the K-band sources observed 
towards the IRDC G192.76+00.10 region. 


3.4. Completeness limit 

The identification of YSOs in a star-forming region using multi¬ 
band data sets is a strong function of the band-pass sensitivities 
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and extinction of the star-forming region. Our primary criterion 
to select YSOs in the present work is the detection in MIPS 24 
pm band. Although the 5cr detection limit of our 24 pm cata¬ 
log is ~ 9.6 mag, the completeness limit of the catalog is ~ 6.7 
mag. We calculated completeness limit of our 24 pm catalog us¬ 
ing artificial stars experiment. In this approach, we insert artifi¬ 
cial stars of different magnitudes to the 24 pm image and then 
perform their detection and phototmetry to recover them. Fig¬ 
ure 5 shows the completeness fraction of artificial stars versus 
corresponding 24 pm magnitudes. This suggests that our 24 pm 
detection is 80% complete up to ~ 6.7 mag. Using SED mod¬ 
els of jRobitaille et al.| ( |2007| ), we roughly estimated that our 24 
pm data is actually more complete towards the high-mass ends. 
The additional YSOs obtained using the K -[3.4] vs. [3.4]-[4.6] 
colour-colour diagram helped us to recover fainter YSOs. The 
WISE catalog is 95% complete at [3.4] = 16.9 mag and [4.6] = 
15.5 mag ( Cutri & et al.||2012) ). We found that 90% of the 3.4 
pm and 4.6 pm sources have a K-band counterpart, thus the ma¬ 
jority of the WISE sources have been detected in our K-band. 
We estimated the completeness limit of our K-band catalog by 
measuring the magnitude at which the cumulative logarithmic 
distribution of sources as a function of magnitude departs from 
a linear slope and begins to turn over. Based on this distribution 
(see Fig. 6), we arrived that our K-band catalog is largely com¬ 
plete down to K ~ 16.5 mag. We tested this approach to our 24 
pm detections, and the resulting completeness limit turns out to 
be ~ 7 mag, comparable to the completeness limit found using 
artificial stars experiment. 


One of the important properties of YSOs is the intrinsic lumi¬ 
nosity. Because of extinction and/or youthfulness of the sources, 
we detected only ~ 50% of the total YSOs in the J-band in com¬ 
parison to the K-band. We thus used K-band magnitudes to es¬ 
timate the luminosities of the YSOs. K-band can be affected by 
the excess emission from the circumstellar material. To account 
for such excess emission, we used average excess emission at 
K-band found in T-Tauri stars. [Meyer et slI.\(1997) found that T- 
Tauri stars typically have a K-band excess between 0.1 and 1.1 
mag, with a median value ~ 0.6 mag. Considering 0.6 mag as the 
excess emission at the K-band for all our YSOs, and using the 
evolutionary model of [Baraffe et al.[ ( [2003] ) for an age of 1 Myr 
(see Sect. 4.1) at 2.5 kpc and at Ay ~ 10 mag, we estimated that 
our K-band completeness level (i.e., K -16.5 mag) corresponds 
to a 0.14 Lq or 0.15 Mq star. Thus, our YSO sample is expected 
to be largely complete above 0.14 L© or 0.15 M©. However, 
for very low-mass stars the adopted excess emission value of 0.6 
mag may not be valid, as in the substellar regime inner part of 
such disks might not emit significant excess emission at K-band. 
Thus, if we consider the emission at K-band only from stellar 
photospheres for low mass YSOs, our completeness limit would 
be - 0.2 L© or 0.2 M©. 


In case of YSOs, it is difficult to derive their true intrin¬ 
sic luminosities without spectroscopic observations. Therefore 
in the absence of short wavelength fluxes or spectroscopic ob¬ 
servations, the above approach appears to be reasonable to have 
approximate luminosity (or mass) of such objects. Since the ef¬ 
fect of excess emission in J-band is minimum, we also estimated 
masses of YSOs using J-band magnitudes. For sources detected 
in both J and K bands, we found that, for majority of our YSO 
the difference in mass estimations is within - 40%. Consider¬ 
ing tentatively that this is the uncertainty associated in the mass 
estimations, the absolute K-band magnitudes of majority of the 
YSO candidates suggest that the region is mainly composed of 
low-mass (< 2M©) YSOs with a mean mass - 0.3 M©. 
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Fig. 7. The spatial distribution of Class I (red stars), Flat-spectrum- 
spectrum (green stars), and Class II (blue stars) YSOs. The continuous 
line running from south-east to north-west marks the highest H 2 column 
density line of the filament along its long axis. 


For the IRDC GI92.76+00.I0 region IRAC four-band data 
sets are unavailable. Thus to quantify the extra percentage 
of Class II and Class I YSOs, we might have detected if 
IRAC observation had been performed, we did analysis of a 
nearby cloud “CI92.75-0.08". The CI92.75-0.08 cloud is lo- 
cated within the S254-S 258 complex and has been studied by 
IChavarria et al.| ( |2008| ) at IRAC bands. We have chosen the 
Cl92.75-0.08 cloud as it is devoid of bright PAH emission like 
the IRDC CI92.76-f00.I0 region, and also at comparable ex¬ 
tinction. To compare, we first selected YSOs of the G192.75- 
0.08 region using the [3.4] vs. [3.4]-[24] diagram as described in 
Sect. 3. 3, and then compared t he statistics with the YSOs identi¬ 
fied by IChavarria et al. (|2008|). W e considered only those YSOs 
from the IChavarria et al.|p0()8| ) catalog whose spectral index is 
greater than -1.6, as done in the present work. This statistical 
comparison yields, if IRAC observation had been performed for 
the IRDC G192.76-1-00.10 region, we possibly would have de¬ 
tected 35% (or 15) more YSOs of Class II, Flat-spectrum, and 
Class I nature. However, it is worth noting that we have already 
added 17 extra YSOs to our MIPS identified YSO sample using 
K - [3.4] vs. [3.4] - [4.6] diagram. This analysis suggests that 
our YSOs statistics in terms of number is comparable, if Class II 
and Class I YSOs had been selected based on only IRAC obser¬ 
vations. 


3.5. Spatial distribution and separation of YSOs 

The spatial distribution of young stars is an useful tool to con¬ 
struct sta£jformation_^enario_ofmobcu]u££loud_compkxes 
(e.g., [Kumar et al. 2QQ7[ "^n et al.|2008t Povich et al.|2009t|De 


March! et al.|2011 

Jose et al. 

20121 Pandey et al. 

2013| 

Mallick 

et al.[[2013 

[Gouliermis et al. 

[20141 [Panwar et a] 

Lp0!4 

[\ [Massi 


et al. 112015 1. In order to understand the star formation scenario 


of the IRDC G192.76-f00.10 region, we have plotted the spa¬ 
tial distribution of its YSO content in Fig. 7. In Fig.7, the con¬ 
tinuous solid line is the highest H 2 column density line of the 
filament along its major axis, and the red, green and blue stars 
are the Class I, Flat-spectrum-spectrum and Class II YSOs, re¬ 
spectively. In order to define the column density line, we first 
created H 2 column density map from the ^^CO column density 
map, and then traced its crests by applying a spatial filtering to 


its intensity distribution. The H 2 column density map was con¬ 
structed from the ^^CO column density map using the relations 
nC^CO)/nC^CO) = 45 SindnC^CO)/n(ll 2 ) = 8 x lO'^ dChavar- 


ria et al.|2008| and references therein). In Fig. 7, one can notice 


that the majority of the YSO candidates are aligned closely with 
the highest column density line of the filament. This close align¬ 
ment strongly suggests that the formation of the YSOs seems 
to be continuing in the dense regions of this filamentary cloud 
along its long axis. We find the mean H 2 column density of the 
1.1 X 10 cm“^, which corresponds to a 


10 mag (using N(H 2 )= 0.94 x 10 Ay cm ^ mag ^ 


highest density line is ^ 
mean Ay 

|Bohlin et al.||I978] ). This is in agreement with the value (i.e, 
Ay ~ 8 mag) expected for filaments above which they are super¬ 
critical and capable of forming stars ( [Andre et al.|20I0| |20I1[ ). 
We find that the projected distances of the identified YSOs from 
the highest column density line range from 0.05 to 1.5 pc, with 
~ 70% falling within 0.4 pc. This narrow separation is a strong 
indication of the fact that the core formation in this filamentary 
dark cloud is not random; their formation has occurred mainly 
along the long axis of the filament. Also these YSOs are proba¬ 
bly formed recently, because they possibly did not have enough 
time to move away from their birth locations. The typical veloc¬ 
ity dispersion seen in rel atively evolved cluste rs and association 
is roughly a few km s“^ ( Madsen et al.|20Q^ . In a recent work, 
[Foster et ^ ( |20I5[ ) showed that the 1-2 Myr old Class II stars of 
the NGC 1333 star forming region have an intrinsic velocity dis¬ 
persion of ~ 1 km s“\ and the average velocity dispersion of the 
dense cores of the region is around 0.5 km/s. Thus, the possibil¬ 
ity that few evolved YSOs might have migrated from their birth 
locations exists, however, in general, protostars are the coldest 
objects, so expected to be embedded within the cores. For exam¬ 
ple, in OMC-3 region, [Thkahashi et al.| ( |2Q13] ) observed eight out 
of sixteen Spitzer sources are associated with the SMA contin¬ 
uum cores (observed with a spatial resolutions ~ 4.5", compara¬ 
ble to MIPS), which are identified as protostars, and six out of 
sixteen are categorized as T-Tauri stars having no counterparts 
in the SMA continuum emission. In the IRDC G192-h00.10 re¬ 
gion, most of the identified YSOs are protostellar in nature. If we 
assume that these protostars have been moved from their birth lo¬ 
cations with a velocity close to 0.5 km/s, their current locations 
(< 0.4 pc) from the highest density line indicate that the forma¬ 
tion of these protostars might have occurred in the last few x 10^ 
yr. 

In young filaments, the distribution of cores is thought to 
represent a possible preferred length-scale of the filament frag- 
mentation (e.g., [Munoz et al. [2007 [[Jackson et al.[2010t[Hacar &| 


[Tafalla|2011t[Miettinen & Qlfner|2013t[Takahashi et al.[[2013[ ). 

Since cores are the precursors to protostars, the distribution and 
separation of young protostars is a good proxy to test the fil¬ 
ament fragmentation theory. To do so, we determined the dis¬ 
tribution of projected separation and projected nearest-neighbor 
(NN) separation among the YSOs as shown in Fig. 8. Although 
the projected separation among the YSOs varies in the range 0.1- 
6.0 pc, the NN separation of the majority of the YSOs shows a 
very narrow range (i.e., 0.1-0.5 pc) with a median ~ 0.19 pc. 
Identification and distribution of cores at longer wavelengths are 
often limited by spatial resolution and sensitivity of instruments. 
So far only a very limited number of millimeter and submil¬ 
limeter observational studies on filaments have been made which 
achieve an angular resolution comparable to, or better than, ex¬ 
isting IR data (i.e., a few arcsec resolution). For example, with 
the high angular resolution observations a separatio n of cores ~ 
0.19 p c ( [Zhang et al.[2009[ beam ~ 1.2"), -0.18 pc dRagan et al. 
2015[ beam ~2"), and ~ 0.25 pc ( [Takahashi et al.[^13| beam 
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Fig. 8. Left: Distribution of the projected distances among the sources 
detected at 24 jum. Right: Distribution of the nearest-neighbor separa¬ 
tion among the YSOs. 


~ 4.5") have been observed in the filamentary clouds IRDC 
G28.34-r0.06, GOl 1.11-0.12 and OMC-3, respectively. These re¬ 
sults seem to be consistent with the median separation of the 
YSOs observed in the IRDC G192.76-h00.10 region. 

Theoretically, it is believed that self-gravitating filaments 
are u nstable to fragmentation and lead to formation of cores 
(e.g., |Nagasaw^|1987[ |Tomisaka]|1995| ). Thus, it is tempting 
to think that the formation of YSOs observed in the IRDC 
G192.76-f00.10 region could be the result of the filament frag¬ 
mentation. In Sect. 4, we therefore discuss whether or not the 
observed distribution and separation of YSOs carry an imprint of 
fragmentation and core formation of the IRDC G192.76-1-00.10 
filamentary cloud. 


3.6. Properties of the YSOs 

In the case of young protostars, it is very difficult to infer their 
stellar and circumstellar properties based on photometric obser¬ 
vations alone. Theoretical models that reproduce the observed 
SEDs provide good representations of the underlying source 
properties. 

To get deep insight into the nature of the detected YSOs, we 
fitted radiative transfer models of |Robitaille et al.| ( |2006|[2n07] ) 
to their observed SEDs. Interpreting SEDs using the radiative 
transfer code is subject of degeneracy and spatially resolved 
multi-wavelength observations can reduce the degeneracy. We 
thus compiled counterparts of the YSOs at NIR (J, H, and K; this 
work), IRAC (3.6, 4.5, 5.8, and 8.0 /tm; |Chavarria et al.|2008|), 
WISE (3.4, 4.6, 12.0, and 22.0 yum; |CutA & et al.|2012| ), MIPS 
(24 pm and 70 pm; this wo rk) and AKARI (65, 90, 140 and 160 
yum; |Yamamura et al.|2009| ) bands, wherever available. We fit the 
models to only those sources for which we have flux values at 
least in five bands in the wavelength range from 1 pm to 24 pm. 
While fitting models to the observed fluxes we adopt the follow¬ 
ing approaches: i) when IRAC 3.6 pm and 4.5 pm fluxes were 
available, the WISE fluxes at 3.4 pm and 4.6 pm were not used; 
ii) for sources not detected at 70 pm or 65 pm, we set the up¬ 
per limits at these bands by assigning the minimum 70 pm or 
65 pm detection flux obtained for the point sources detected in 
the region. The 70 pm flux upper limit is assigned to only those 
sources that are not detected in the 70 pm observed area. Sources 
that are located outside the 70 pm observed area and not detected 
in 65 pm band, we assigned the minimum 65 pm flux obtained 
for the point sources observed in the region from the AKARI 
survey; iii) we scaled the SED models to the distance of the fila¬ 
ment (i.e., 2.5 + 0.2 kpc; [Chavarria et al?] ( |2008| )) and allowed a 
maximum Ay value determined by tracing back the YSOs’ cur¬ 
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rent location on J vs. J-H or K vs. H-K or K vs. K-[3.4] diagram 
to the intrinsic dwarf locus along the reddening vector (see e.g., 
ISamal et al.|20T0l i. 


Eigure 9 shows the SEDs of 38 sources that satisfy our five 
data point flux criteria. Due to lack of optical, far-infrared, and 
millimeter data points, it is quite apparent that the SED models 
show high degree of degeneracy; nonetheless, the SEDs clearly 
indicate that the majority possess IR-excess emission, possibly 
emission from circumstellar disk and envelope. While looking 
at Eig. 9, we found that the observed SEDs of four sources (IDs. 
4, 18, 19 and 29 in Eig. 9) are fitted well by reddened stellar 
photosphere models. We find that they correspond to those four 
sources that have been rejected as YSO candidates in Section 
3.2 on the basis of their location on the [3.4] vs. [3.4]-[24] di¬ 
agram (i.e., sources with [3.4]-[24] colour < 1 mag). Thus, the 
SED models confirm our previous results that the bright 3.4 pm 
sources found in the field-star/MS zone of [3.4] vs. [3.4]-[24] 
diagram are indeed disk-less stars. 


It is not possible to characterize all the SED parameters from 
the models du e to limited observation al data points. However, 
as discussed in [Robitaille et al.| ( |2007| ), some of the parameters 
can still be constrained depending on the available fluxes. Eor 
example, in the present study, for majority of the sources, the 
SED models between 1 pm to 70 pm represent fairly well the 
data points, hence, the disk parameters are expected to be better 
constrained. However, it is worth noting that precise determi¬ 
nation of disk parameters using SED models, as demonstrated 
by jSpezzi et al.j ( |2013| ), requires data from optical to millime¬ 
ter bands, as well as good knowledge of the p hysical parame- 
ters of the central star. Similarly, as discussed in [Robitaille et al.j 
( [2QQ6| ), in the case of young protostars, the disk is expected to 
be deeply embedded inside the envelope and the relative contri¬ 
butions of the disk and envelope to the SED are difficult to dis¬ 
entangle. Because of the above reasons, in the present case, even 
though the uncertainty in the disk parameters is expected to be 
high, it should nonetheless give a good proxy of the nature of the 
sources. In Table 2, the disk parameters (disk mass (Mdisk) and 
disk accretion rate (Mdisk)) are tabulated. Since our SED mod¬ 
els are highly degenerate, the best-fit model is unlikely to give 
an unique solution. So, the tabulated values in Table 2 are the 
weighted mean and standard deviation of the physical parame¬ 
ters obtained from the models that satisfy - Ymin - ^^data 
weighted by of each model, where is the goodness- 

of-fit parameter for the best-fit model and Adata is the n umber 
of input observational data points (see [Samal et al.|2QT^ . Erom 
Table 2, we found that the disk masses and disk accretion rates 
of ~ 80% YSOs are in the range 0.001 - 0.021 Mq and 0.01 - 
0.95 X 10“^ Mo/yr, respectively. We note that in addition to the 
limitations outlined above, the absolute uncertainties associated 
with the disk p arameters are possibly in a range of 1-3 orders of 
magnitude (see [Robitaille et al.|2007[ ). Therefore, we stress that 
the derived disk parameters must be considered as representative 
values, and should be treated with caution. 


Due to several limitations, it is not possible to comment on 
individual objects, but if we take the SED results for statistical 
purpose, the derived disk accretion rates suggest that the YSOs 
of the IRDC G192.76-f00.10 region are mainly low-mass in na¬ 
ture (e.g., [Dahm|2008[ ), consistent with the nature of the sources 
derived from the photometric data. 
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4. Discussions and Conciusions 

4.1. Evolutionary status 

Protostar fraction (i.e., number of protostars (Class I + Flat- 
spectrum) of the total number YSO population) in young clus¬ 
ters is a good tracer of age. For example, the protostar fraction 
derived from the data involving 2.2-24 jim is 14% in the IC 
348 cluster of a ge ~ 2-3 Myr, ~ 16% in Chamaeleon II ([Al¬ 
cala et al.|2008| ) of age ~ 2-3 Myr ( Sciortino|2007[|Spezzi et al. 


per unit length) valu e exceeds its cri tical equilibrium mass value 
Merit = '^c'^sIG (e.g., Ostriker 19641, where c^ is the sound speed 


2008|), and ~ 36% in the NGC 1333 cluster of age ~ 1-2 Myr 
( |j0rgensen et al.|[2Q06| ). This suggests an evolutionary differ¬ 


ence, with NGC 1333 being young consists of more younger 
population of YSOs compared to IC 348 and Chamaeleon II. 
In the IRDC G192.76-f00.10 region, we have identified a to¬ 
tal of 62 YSOs, of which 19 are Class II YSOs and 43 are 
Class I plus Flat-spectrum YSOs. This indicates that this dark 
cloud is unusually rich in protostars with fraction ~ 70 %. This 
high fraction strongly suggests that the IRDC G192.76-f00.10 
region is too young for most YSOs to have reached the Class II 
stage. In young clusters, the error in the Class ratio in general 
is dominated by the detection of the actual number of Class I, 
Flat-spectrum and Class II objects. If we consider all the YSOs 
above our completeness limit, we find that the protostar frac¬ 
tion is still high (i.e., ~ 55 %), indicating that the region indeed 
contains a high percentage of protostars. A comparison of the 
protostar fraction of the IRDC G192.76-f00.10 region with the 
IC 348, NGC 1333, and Chamaeleon II suggests that the IRDC 
G192.76-1-00.10 is indeed young, possibly younger than 1 Myr. 

Another possible way to assign age to a star-forming region 
is to use the lifetime of different phase of YSOs associated in 
the region. [Evans et al.| ( |2009| ) derived lifetime of different YSO 
phases from a large sample of YSOs collected from five nearby 
molecular clouds, assuming that in these clouds star formation 
has proceeded at a constant rate over 2 Myr period of time and all 
the prestellar cores have evolved into Class O/I YSOs, then into 
Flat-spectrum YSOs and then to Class II YSOs. Although the de¬ 
rived lifetimes are still subject to large uncertainties (see [Evan^ 
jet al. 12009) and discussion therein), the above work involving a 
large set of YSOs suggests that in general a lifetime of Class II, 
Class I, and Elat-spectrum stages of a YSO is roughly ~ 2, ~ 
0.44, and ~ 0.35 Myr, respectively. In the IRDC G192.76-1-00.10 
region, we find that most of the identified YSOs are Class I, Elat- 
spectrum, and Class II in nature. If we use the above lifetimes, 
then the YSOs class statistics suggest a mean age of ~ 1 Myr for 
IRDC G192.76-f0.10. 


4.2. Possible fragmentation process of the IRDC region 

The S254-S258 complex harbors a 20 pc long filament. Nu¬ 
merical simulations suggest that there are various processes by 
which filaments can be formed. Simulation shows that the frag¬ 
mentation of clouds into sheets and filaments are the natural 
consequence of supersonic turbulence present in the inter-stellar 
medium (jKlessen & Burkert|2000t[Padoan et al.|2001[|Ostriker| 


et al.|200T Bate|2009 ). The driving sources for large-scale tur- 

bulence could be flows of atomic gas (e.g., [Hennebelle et al.j 
2008[[BSierjee et al.|2009|) or waves from supernova explosions 


and superbubbles (e.g., |Matzner|2002[|Dale & Bonnell|2011|) o r 

collisions between molecular clouds (e.g.. Tasker & Tan|2009) . 
Once filamentary structures formed, theoretical models suggest 
that they are subject to fragmentation (e.g. , |Larson| 1985 ). 

Eor an infinite and isothermal filament, the filament is unsta¬ 
ble to axisymmetric perturbation, if its line mass Mune (i.e., mass 


of the medium and G is the gravitational constant. The critical 
line mass only depends on the gas temperature. The temperature 
of the IRDC G192.76-1-00.10 cloud is ~ 14 K (discussed below), 
which corresponds to Merit ~ 25 Mq pc“^. 

The projected size of the IRDC G192.76-r00.10 cloud along 
its long axis is about 5.7 pc. To estimate Mune, we first estimated 
the radial profile of the filament from the H 2 column density 
map. To do so, at several positions along the filamentary struc¬ 
tures, we extracted perpendicular column density profiles and 
established the radial extent of each profile using Gaussian fit¬ 
ting similar to other works (e.g., [Arzoumanian et al.|20TT[ [Smith 


et al. 


2Q14[ ). We then estimated Mune by dividing the length of 


the filament to the mass estimated over the mean radial extent 
of the entire filament, which turns out to be ~ 120 Mq pc“^ 
The observed line mass is ~ 5 times larger than the critical line 
mass, indicating that the filament is supercritical, hence suscep¬ 
tible to fragmentation. Supercritical filaments are believed to be 
globally unstable to radial gravitational collapse and fragment 
into prestellar clumps or cores along their major axis. This is 
supported by the roughly linear sequence of protostars observed 
along the long axis of the IRDC G192.76-h00.10 region. Recent 
Herschel observations have also shown that supercritical fila¬ 
ments usually harbor several prestellar cores and Class 0/Class I 
protostars along their length (e.g., [Andre et al.|201Q[ ). 

Numerical simulation also suggests that gravitational frag¬ 
mentation is possibly the dominant mechanism of formation of 
cores and stars in the dense regions of the filaments (e.g., [Klessen[ 
et al.|2004| ). To further validate the above hypothesis, we com¬ 
pare our results with the simple mo del of gravitationa l fragmen- 
tation, given by [Hartmann[ ( 2002| ). In the model of [Hartmann[ 
( 2002| ), the fragmentation length (T^) of a cylinder due to gravi¬ 
tational instability is given by: 


Ac = 3.9Ac]l{GJ.) = 


( 1 ) 


where. Ay is the visual extinction corresponding to the surface 
density (E) through the center of the filament and Tio is the 
gas temperature in units of 10 K. The corresponding collapse 
timescale is given by: 


T ~ 3.7r;/'A^'Myr. 


( 2 ) 


Herschel estimates of dust temperatures indicate that the 
temperature of the quiescent and star-forming filaments are in 
the range of 10 K to 15 K ( [Arzoumanian et al.|2011| ). Observa¬ 
tions of the ammonia (1,1) and (2, 2) inversion transitions sug¬ 
gest that the kinetic temperature in the dense part of the IRDC 


G192.76-1-00.10 filament is in the range of 13 to 15 K (Dunham 
[et al.|2010| . Considering 14 K as the temperature of the filament 
and 10 mag as the mean visual extinction, the above equation 
predicts Xe ~ 0.21 pc and r ~ 0.45 Myr. Eor randomly oriented 
filament, the line of sight median inclination angle could be ~ 
60° ( [Genzel & Stutzki|1989| ). If this is the case, then this would 
decrease Xe to a value 0.19 pc. Most of the YSOs in the filament 
are protostellar in nature. The lifetime of embedded protostars 
is around 0.5 Myr. The predicted fragmentation length and col¬ 
lapse time for the region are comparable to the observed median 
NN separation (~ 0.2 pc) and age of majority of the YSOs; thus 
the model predictions are broadly in agreement with the obser¬ 
vations. 
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4.2.1. Limitations 


The scenario pres ented above is rough l y con sistent with the re¬ 
sults discussed in [Ballesteros-Parede^ ( |2004| ) and with the sce¬ 
nario preferred for the YSOs formation in the Taurus ( |Hartmann| 
|2002| ). However, there are number of issues yet to be addressed. 
For example, the exact role of turbulence prior to fragmenta¬ 
tion is not known to us; its presence would increase the effec¬ 
tive sound speed, thus would increase the critical mass per unit 
length. However, it is suggested that pre-stellar cloud is most- 
likely stabilized against global collapse by interstellar turbu¬ 
lence, but this support is quickly removed on small spatial scales 
(e.g., [Elmegreen et al.||200Q| ). Recent observations suggest that 
the core formation in filamentary cloud is possibly a two-step 
process, in which first, the subsonic, velocity-coherent filaments 
condense out of the more turbulent ambient cloud. Then, the 
cores fragment quasi-statically, where turbulence seems to play 
a little or no role in the formation of the individual cores (e.g., 
Ballesteros-Paredes||2004t [Andre et al.||2011t |Hacar & Tafalla] 


201 Ij ). Based on ammonia observations, [Dunham et al.n2010 ) 
observed non-thermal velocity dispersion of -0.43 km s“^ in the 
central region of IRDC G192.76-h00.10. However, non-thermal 
dispersion in clouds could be of various origin. For example, 
non-thermal motions can be driven by the accretion flows to 
the filament potential (e.g., [Klessen & Hennebelle|2010[ ), if the 
IRDC G192.76-r00.10 filament has undergone dynamical evo¬ 
lution after fragmenaion. Similarly, protostellar feedback can 
also drive supersonic turbulence at sub-parsec size scales (e.g., 
[Duarte-Cabral et aT][2012[ ). If the non-thermal motion (corre¬ 
sponding to the velocity dispersion) is to be considered in addi¬ 
tion to the thermal component to support the filament, this would 
then increase the effec tive sound speed of the filament (e.g., see 
MacLaren et al.|19^ . As a consequence of this, the Merit of the 
filament would increase to ~ 75 M© pc“\ which is still lower 
than than the observed Mune, implying that the total amount of 
support available is still insufficient to keep the filament in equi¬ 
librium, hence the filament should be radially contracting. 

Other physical processes that might affect the model parame¬ 
ters discussed above, are helical magnetic field ( [Fiege & Pudritz[ 
2000[ ) and star formation due to dyna mical effects such as flows 


from surrounding sub-filaments (e.g., [Smith et al.|20I4[[Schn^ 
der et al.|20I0[ ). Helical magnetic fields are believed to decrease 


the cr itical length of fragmentation of a cylinder ( [Fiege & Pudritz[ 


20001, whose effect in the present case is difficult to assess. How¬ 


ever, in case of star formation due to converging flows, higher 
concentration of dense gas is expected at the junction point of 
the converging sub-filaments and the main filament, which tend 
to increase star formation activity. For example, in the DR21 
region several low-density striations or sub-filaments were ob¬ 
served perpendicular to the main filament and apparently feeding 
matter to t he main filament (e.g., [Kumar et al.|2007[ [Schneider] 


[et al.|2010 [Hennemann et al.|2012[ ). Considering that the IRDC 

G192.76-r00.10 region is slightly structured at its eastern side, 
the possibility that few stars might have formed by other dynam¬ 
ical processes can not be ignored. However, the morphology of 
the IRDC G192.76-h00.10 area is largely linear in CO and we do 
not see strong morphological signature of such structures (e.g., 
thin long sub-filaments) that are radially attached to the IRDC 
G192.76-h00.10 region, though this cannot be warranted without 
high resolution and sensitive observations. Nonetheless, to limit 
the impact of dynamical processes on our estimations, we cal¬ 
culated again the model parameters, without accounting for the 
group of stars seen around (r2ooo=93.50 & ^ 2000 = +17.87 (i.e., 
the group of stars seen distinctly away from the filament’s long 
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axis), as well as a few stars located far away from the highest 
density line (see Fig. 7). Considering tentatively that these are 
stars that might have formed by some other processes, we found 
that the model parameters are changed by 20% only. 

We also stress that even though smooth cylindrical models 
are useful to understand star formation processes in filaments, 
they can only represent a first order approximation because real 
clouds are likely to have much more density inhomogeneities 
than what we have assumed here. In addition to above, the initial 
cloud configuration also plays a decisive role where stars would 
form. For example, if the original parental cloud is not infinitely 
long or spherical, gravitational edge focusing results in enhanced 
concentrations of mass at one end of the filament, where local 
collapse would proceed faster than collapse of the entire fila¬ 
ment (e.g., [Burkert & Hartmann|2004[ ). We do not observe high 
concentration of stars or massive condensations at any end of the 
IRDC G192.76+00.10 region, thus the gravitational edge focus¬ 
ing effect is not preferably happening here. 

In summary, the caveats outlined above prevent us to make 
any definite statement on the YSO formation; however, if we 
consider the role of the helical magnetic field, and star forma¬ 
tion by other dynamical processes are minimum, then we find 
a remarkable reconciliation of the observed properties with the 
model predictions, suggesting that gravitational fragmentation 
is probably the dominant cause of core formation in the IRDC 
G192.76+00.10 filamentary cloud. Confirming and refining the 
scenario require high angular resolution mid-infrared, dust po- 
larimetric and millimeter line observations to identify other pos¬ 
sible low-mass YSOs, to constrain roles of magnetic field, and 
influence of gravity, respectively in the star formation processes 
of the region. 

4.3. Star formation efficiency and rate 

In this section, we discuss the star formation efficiency (SFE) 
and star formation rate (SER) of the IRDC G192.76+00.10 re¬ 
gion based on the YSOs identified in the present work. SEE and 
SER are the fundamental physical parameters that are essential 
for the understanding of evolution of star-forming regions and 
galaxies. 

SEE is defined as the ratio of the total stellar mass to the total 
mass of stars and gas. We estimate the total gaseous mass (M^^^) 
associated to the IRDC region -1100 M©. To estimate the 
we used the ^^CO column density map and integrated the column 
density four times above local background, i.e., column density 
value > 1.1 X 10^^ cm“^ over 11 pc^ area. We then converted the 
total column density to the total H 2 column density as dis¬ 
cussed in Sect. 3.5, and used the following equation to estimate 
the mass: 


M = yumHApixSH2 (3) 

where, p is the mean molecular weight, mn is the mass of the hy¬ 
drogen atom, S H 2 is the summed H 2 column density, and Apix 
is the area of a pixel in cm“^ at the distance of the region. The 
determination of the stellar mass in accreting stars is not trivial, 
as emission due to circumstellar disk surrounding the star and 
accreting material from the protoplanetary disk onto the central 
star affects the observed spectrum. Broad-band spectroscopic ob¬ 
servation would reveal more accurate mass of the YSOs (e.g., 
[Manara et al.|2013] ); however, if we assume that each source of 
the IRDC G192.76+00.10 region has a mass of 0.5 M©(close to 
the mean stellar mass of the region; see Sect. 3.4), consistent 
with the characteristic mass from the studies of [Chabri^ ( [2003[ ) 



















































M. R. Samal et al.: A cluster in the process of formation 


and |Krou^ ( |2001 1 ) type of initial mass function (IMF), then the 
total mass of all the YSOs (Mysos) is ~ SOM©. Using the 
and we estimated the SFE ~ 0.03 or 3% in the IRDC 

G192.76+00.10 region. . 

SFR describes the rate at which the gas in a cloud is turning 
into stars. Assuming 1 Myr as the duration of star formation (see 
Sect. 4.1) and using the derived and we estimated 

the SFR (= xSFE/tsf, where tgf is the star formation time 
scale) of the IRDC G192.76+00.10 region ~ 30 M© Myr“^ The 
projected area of the IRDC G192.76+00.10 region over which 
we estimated the cloud mass is ~ 11 pc^. If we normalized the 
derived SFR by the cloud area, this leads to SFR per unit sur¬ 
face area ~ 3 M© Myr"^ pc“^. The SFE and SFR per unit sur¬ 
face area of the IRDC G192.76+00.10 region are comparable 
to the SFE (i.e., 3-6%) and SFR per uni t surface area (i.e., 0.6- 
3.2 M© Myr“^ pc“^) values reported by Evans et al. (2009) for 
nearby molecular cloud complexes. 

Since our YSOs sample is complete down to ~ 0.15 M©, 
thus we may be missing a population of deeply embedded YSOs 
of masses below this completeness level. If we assume that the 
cloud has already formed stars down to the hydrogen burning 
limit (i.e., ~ 0.08 M©), this would not drastically alter the SFE 
of the complex, because most of the observations in our Galaxy 
are consistent with an IME that declines below 0.1 M©(e.g.,|Lada 


2005 [Oliveira et al.||2009t [Andersen et al.||2008t [Neichel et al. 


2015). Thus, we do not expect a large population of Class II and 
Class I YSOs below our completeness limit. The shape and uni¬ 
versality of the IME at the s ub-stellar regime is still under debate 
though, however, assuming Chabrier[ ( [2003 [ ) type of mass func¬ 
tion, we may probably be missing only about 20% of the total 
number of YSOs. This may increase the SEE of the region by 
only about 1%, assuming each YSO has a mass 0.5 M©. 


4 . 4 . Emerging young cluster 


The detection and identification of clumps and cores in star¬ 
forming complexes using long wavelength observations depend 
actually on the resolution and sensitivity of surveys. Eor exam¬ 
ple, the dust clump “PI” of the G28.34+0.06 region detected 
with the IRAM 30 m telescope and JCMT at resolutions from 
ll"to 15''( [Rathbome et al.|2006[[Carey et al.|1998 ) are resolved 
by the Submillimeter Array (SMA) into five cores at 1".2 resolu¬ 
tion ( [Zhang et al.|2009[ ). So far, no high resolution observations 
have been conducted to identify clumps and cores of the IRDC 
G192.76+00.10 region. In this work, we have detected protostel- 
lar sources of mass down to 0.15 M©, so millimeter continuum 
observations of mass sensitivities better than 0.45 M© (assuming 
SEE of cores is ~ 30%; [Alves et al.[ ( [2007[ )) and spatial resolu¬ 
tion better than Spitzer observations would be very helpful for 
the detection of resolved cores in the region. Nonetheless, using 
1.1 millimeter Bolocam Ga lactic Plane Survey (BGPS; beam ~ 
33"), Dunham et al. ( 2010 ) observed an elongated clump (their 
ID 50) at the center area of IRDC G192.76+00.10. Prom their 
observation, we estimated the virial parameter ( [Bertoldi & Mc-[ 
[Kee[[1992| ) of the clump -0.3. Virial parameter represents the 
dynamical stability of a core or clump. A clump is assumed to 
be gravitationally bound if its virial parameter value is < 1, un¬ 
bound otherwise. The virial paramete value of the BGPS clump 
suggests the clump is susceptible to gravitational collapse, and 
can lead to star formation. 


The filament in which the BGPS clump is embedded, is 
in super-critical (discussed in Sect. 4.2.1) stage, implying that 
the filamentary cloud may be contracting. Moreover, observa¬ 


tions and simulations of filaments show accretion fiows along 
the filaments onto the clumps and cores located at the bottom 
of their potential ([Balsara et al.[[2001t [Kirk et ar][2Q13t [Tack-[ 
[enberg et al.[[2014[ ), and some clumps and cores possibly ac¬ 
crete mass faster than others due to their higher gravitational 
potential (e.g., [Tackenberg et al.[[2014[ ). Observations suggest 
that star-forming clouds dynamically evolve on time-scales of 
a few Myr. If it is to be believed that cluster formation mainly 
occurs in filamentary clouds, then the lifetime of the typical fila¬ 
ment should be a few Myr, since clusters with age greater than 5 
Myr are found to be seldom associated with molecular gas (e.g., 
[Leisawitz et al.[[1989[ ). In Sect. 4.1, we have estimated that the 
age of the IRDC G192.76+00.10 region is around 1 Myr. If we 
put in context all our results, they suggest that the YSOs of the 
IRDC G192.76+00.10 region are young, they are embedded in a 
gas reservoir of mass - 1100 M©, and the dark cloud is forming 
stars at a high rate of - 30 M© Myr“^ and is part of a large fila¬ 
ment environment where individual clumps and cores can grow 
in mass by accreting matter from the parental filament to their 
potentials. These evidences suggest that if the dark cloud will 
evolve dynamically for a few Myr, then the possibility that it 
may emerge to a richer cluster exists. 


4 . 5 . Overall picture of star formation In the complex 

The 20 pc long filament of the S254-S258 complex consists of 
six evolved H ii regions. Among these, four H ii regions (S254, 
S255, S256, and S257) are located at the center of the complex 
at a mean separation of 2.6 pc (see Pig. 1). The exciting stars of 
S254, S255 and S257 are isolated sources, without any co-spatial 
clustering of low-mass stars around them. The mean dynamical 
age of these evolved H ii regions is - 2.5 Myr (see [Chavarria| 
[et al.[|2008) ). It seems that the evolved massive stars (exciting 
stars of evolved H ii regions) in this complex must have formed 
about 2.5 Myr ago. 

In this work, we characterized the protostellar population of 
the IRDC G192.76+00.10 region, which is located farther away 
from the evolved H ii regions of the complex. Although the low- 
mass YSOs in the IRDC G192.76+00.10 region are at different 
evolutionary stages, the majority (~ 70%) of them are Class I and 
Plat-spectrum in nature. The typical lifetime of Class I and Plat- 
spectrum objects is around 0.4 Myr ( [Evans et al.|2009[ ). There¬ 
fore, it is very likely that the majority of the YSOs in the IRDC 
G192.76+00.10 region were formed simultaneously or over a 
narrow range of time, possibly less than 0.5 Myr of time. 

The evolutionary status of the hig-mass OB stars (age -2.5 
Myr) and low-mass protostars of the IRDC G192.76+00.10 re¬ 
gion (age - 0.5 Myr) suggests existence of multi-generation star 
formation in the complex. The evolved H ii regions (e.g., S255, 
and S257) appear to be projected on the long filament. If all the 
optically visible evolved H ii regions and the low-mass proto¬ 
stars of the IRDC G192.76+00.10 region are part of the same 
long (-20 pc) filament, it is then plausible to think that multiple 
generations of stars may be residing in the whole filament. If this 
is the case, then the whole filament probably has not fragmented 
as a single entity. Hierarchical fragmentation has also been ob¬ 
served in several filamentary clouds (e.g., [Hacar&Tafalla[2011t 
Takahashi et al. 2013[). Eor example, in Orion filamentary cloud, 
Takahashi et al. ( [2013[ ) found different fragmentation scales for 
large-scale clumps (9-10 pc), small-scale clumps (- 2.5 pc), and 
dense cores (0.15-0.55 pc), and interpreted this result as the hi¬ 
erarchical fragmentation of Orion A filament. Thus, the possi¬ 
bility that the long filament of the S254-S258 complex has un¬ 
dergone hierarchical fragmentation can not be ruled out. Hier- 
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archical fragmentation in clouds could be of various origins. 
For example, it could be due to combined effect of turbulence, 
cloud temperature, magnetic field, and ionization from the mas- 
si ve stars as discussed in|Takahashi et al.| ( |2Q13| ) or as discussed 
in |Vazquez-Semadeni et al. (2009), where low-mass and high- 
mass star forming regions can arise from hierarchical gravita¬ 
tional collapse of large cloud, where the former regions arise 
from collapse of small-scale fragments, while the latter may ap¬ 
pear due to large-scale collapse. 


4.6. Concluding remarks 

In this paper, we have presented a dark cloud (mainly consists 
of low-mass stars) based on MIPS, WISE and JHK photome¬ 
try. It harbors 62 YSOs, 49 of which contain 24 pm sources, 
indicating that the cluster is young. Interestingly, the dark cloud 
is found in the filament of an OB complex, with members dis¬ 
tributed in an aligned fashion along the long axis of the fila¬ 
ment and at a NN separation ~ 0.2 pc. We find that the pro¬ 
tostars to Class II ratio in the dark cloud is high (~ 2). When 
compared to other star-forming regions, it indicates that the re¬ 
gion is young, possibly younger than 1 Myr. From our results, 
we also find evidences favoring gravitational fragmentation as 
the dominant cause of YSOs formation in this dark cloud, but a 
conclusive answer would require a further large-scale study of 
the filament concerning magnetic field, turbulence, and gravity. 
The lack of deep infrared observations prevent us making a defi¬ 
nite conclusion on the properties of the cloud; however, from the 
present observations, we argue that the cloud is currently form¬ 
ing a cluster at an efficiency of ~ 3% and a rate of ~ 30 Mq 
Myr“^. The cluster is embedded in a supercritical filamentary 
dark cloud of mass -1100 Mq and is part of a large filament 
geometry. We hypothesized that it may become a richer cluster 
with time. Due to the lack of massive protostars, the possibility 
of radiative heating and photo-evaporation of the dense gas in 
the studied filament region is low, hence the dark cloud provides 
an unique opportunity with high sensitive observations to gain a 
deeper understanding of cluster formation and evolution within 
the filamentary environment. 
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Table 2. Inferred Physical Parameters from SED Fits to YSOs 


ID 

JJdisk 

(Mo) 

Mdisk 

(lO-*^ Mo/yr) 

01 

0.006 + 0.01 

0.10 + 0.22 

02 

0.007 + 0.02 

0.08 + 0.27 

03 

0.013 + 0.02 

0.64+1.49 

05 

0.006 + 0.01 

0.24 + 0.66 

06 

0.007 + 0.01 

0.31+0.78 

07 

0.009 + 0.02 

0.23 + 0.72 

08 

0.007 + 0.01 

0.09 + 0.60 

09 

0.012 + 0.02 

0.95 + 1.43 

10 

0.025 + 0.04 

3.10 + 3.38 

11 

0.003 + 0.01 

0.01 + 0.04 

12 

0.030 + 0.03 

5.79+12.35 

14 

0.006 + 0.02 

0.14 + 0.47 

15 

0.005 + 0.02 

0.47 + 2.14 

16 

0.043 + 0.04 

3.88 + 5.48 

17 

0.007 + 0.01 

0.21 + 0.47 

20 

0.006 + 0.01 

0.34+1.37 

21 

0.041 + 0.03 

4.63 + 3.58 

22 

0.005 + 0.01 

0.65 + 1.03 

23 

0.002 + 0.01 

0.02 + 0.09 

24 

0.002 + 0.01 

0.02 + 0.09 

25 

0.078 + 0.09 

0.48 + 0.79 

26 

0.014 + 0.03 

2.96 + 3.99 

28 

0.016 + 0.02 

1.46 + 2.38 

31 

0.005 + 0.01 

0.11+0.29 

32 

0.006 + 0.01 

0.40 + 0.97 

33 

0.015 + 0.02 

1.89 + 2.67 

34 

0.005 + 0.01 

0.14 + 0.26 

36 

0.015 + 0.02 

0.37 + 0.85 

42 

0.021 + 0.02 

1.42 + 2.37 

44 

0.015 + 0.02 

0.32 + 0.73 

45 

0.018 + 0.04 

0.48 + 1.09 

47 

0.029 + 0.03 

8.88 + 14.02 

48 

0.005 + 0.01 

0.18 + 0.50 

49 

0.001 + 0.02 

0.13 + 0.02 
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Fig. 9. SEDs of the YSOs discussed in the text. The black line shows the best fit model, and the grey lines show subsequent models that satisfy 
- 2^data Criteria. The dashed line shows the stellar photosphere corresponding to the central source of the best fitting model. The circles 
denote the observed fiux values. The IDs correspond to the source given in Table. 2. 
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